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ABSTRACT 



We compute the strong coupling constant gp*p n , with P and P* respec- 
tively pseudoscalar and vector heavy mesons by using the QCD sum rules 
approach. Our computation is based on the evaluation of the time ordered 
product of currents between the vacuum and the soft pion state. The so-called 
parasitic terms are taken into account and give a contribution to the sum rule 
of the same order of the lowest lying state, while higher dimension non per- 
turbative terms have small numerical effects. The infinite heavy quark mass 
limit is also examined. 



1 Introduction 



In this letter we shall address the problem of the determination of the strong coupling 
constant gp*p n , where P and P* are respectively the lowest pseudoscalar and vector heavy 
meson made up of a light q antiquark and a heavy quark Q. This coupling is defined by 
the strong amplitude 

< 7r-(g) P°(q 2 )\P*-{ qi ,e) > = g P * Pn e» ■ q, (1.1) 

and is of physical interest for a number of reasons. First, for P = D and P* = D*, the 
decay amplitude Ql.lQ describes a process which is observed experimentally. Present data 
provide the upper bound: 

9d*Dw < 13.8 (1.2) 

from T tot (£>*+) < 131 KeV gj and BR(D* + -> D°tt + ) = 30.8 ± 0.4 ±0.8 §. Second, 
it is commonly believed that the form factor F\(q 2 ) describing the semileptonic decay 
B — ► 7r£z/ is dominated, for q 2 ~ g^ aa . = {m B — m„-) 2 , by the 5* pole; in this case Fi(q 2 
would be given by the formula H. M: 



lmax J 



2 , 111" VIS' Hr 



* T m , + (m a .-m B ) (L3) 

so that knowledge of the strong coupling gB*B-w would shed light on this decay process 
(fs* is the leptonic decay constant of B*). Finally in the chiral effective field theories for 
heavy mesons that have recently received attention by a number of authors ||, 0, 0, 
||, the above mentioned coupling is related to the interaction lagrangian Q 

C mt = ~ g Tr[H llll5 {t}d^ - ^)H] (1.4) 



by the formula 

The coupling constant g appears in a number of calculations that make use of the chiral 



gp*Pn = —^-g ■ (1.5) 



effective theories in the infinite heavy quark mass limit [§,0, [ 10 1 , [11 



The computation we shall describe in the next sections is based on the QCD sum 
rules approach 0. Some other QCD sum rules calculations of the strong coupling constant 
g P * Pn can be found in the literature [0, fL4|. The authors of ref.[[H| compute go*D-n by 



considering two-point functions in an external pion field, whereas the authors of ref.fH} 
compute the coupling constant g of eq. ( P-.4|) . We differ from these calculations in several 
aspects. First of all we compute gD*Dn, gB*B-w and g, which allows us to discuss l/mq 



corrections to the asymptotic itlq — > oo limit. Second we use the method of ref.||15|| both 
in the evaluation of the correlators (the time ordered product of currents is taken between 
the vacuum and the pion state) and in the treatment of the so called parasitic terms, 
arising from non-diagonal transitions Q they are present if, as in our case, one takes the 

1 We employ the notations of where H = (l+ft)/2(P* 7 M — P575), (v is the heavy meson velocity, P* 
and P5 are annihilation operators of the heavy mesons P* and P respectively) and £ 2 = S = expiM/f, 
with M the 3x3 matrix of the Nambu-Goldstone pseudoscalar meson octet and / = f n in the chiral 
limit {f n = 132 MeV). 

2 For a review of QCD sum rules as applied to heavy masses see |^2| . 

3 See also @ " 
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soft pion limit q — > 0, which implies q\ = q2 and one has to take the single and not the 
double Borel transform. These differences alter the numerical result, as we shall discuss 
below; other differences with refs. [13|, [Q, arising from higher order non perturbative 



terms that we have included, have small numerical effects in the final results. 



2 QCD sum rule calculation of gD*Dn and gB*Bn 



For definiteness we shall consider the off-shell process 

B*-(e, qi )^B°(q 2 )+7r-(q), 
the case of D* — > Dtt being completely analogous. We consider the correlator 



A,(P,q) =i dx< 7r-(g)|T(^(x)j 6 (0)|0 > 



-iqix 



Aq^ + BP^ 



'XI) 



(2.2) 



where = wy^b, js = ib^d, P = qi + q2 and A, B are scalar functions of qf, q 2 , q 2 . 

Both A and B satisfy dispersion relations and are computed, according to the QCD 
sum rules method, in two ways: either by saturating the dispersion relation by physical 
hadronic states or by means of the operator product expansion (OPE). We shall focus 
our attention on the invariant function A that we compute in the soft pion limit (q — > 0) 
and for large Euclidean momenta (q\ = q\ — > — oo). 

Performing the OPE on A we get the following result 

A = A® + + A^ + A® + A^ + A® (2.3) 



with 



A {1 
A^ 
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2m 2 



(ql - m 2 b ) 2 (ql - m 2 b f 



(2.4) 



In eqs.( [2.4| ) < uu > is the quark condensate (< uu >= — (2A0MeV) 3 ), m and mi are 
defined by the equations 

< ug s cr ■ Gu >= ml <uu> (2.5) 



and 



< Tr(q)\uD 2 j^j 5 d\0 >= -if^m 2 ^ 
2 



(2.6) 



and their numerical values are: wig 



0.8 GeV 2 , m\ = 0.2 GeV 2 [|TJ, |T7J. The origin of 
the different terms in (|2.3| ) is as follows. is the leading term in the short distance 
expansion; and A^ arise from the expansion of V^(x) at the first, second and 

third order in powers of x; and A® arise from the expansion of the heavy quark 
propagator at the second order and from the zeroth and first term in the expansion of 
Vfj,(x) respectively. We stress that we have considered all the operators with dimension 
D < 5 in the OPE of the currents appearing in ( |2.2|) ; we also note that, with the obvious 
change b — > c, eq.(j2.4j) also applies to the amplitude for charm decay D* — ► Dir. 



Let us now evaluate the hadronic side of the sum rule, which can be obtained by 
writing for A(0, qf, q%) the dispersion relation 

A(0,q 2 ,ql) — — J dsds 2\T~> 2\ + subtractions . (2.7) 

7T J (s — qi)(s — q 2 ) 

We divide the integration region into three parts. The first region (I) is the square given 
by ml < s < s , mf < s' < s ] for s small enough, (I) contains only the B and B* poles, 
whose contribution is 



Ai(0,ql,q 2 2 ) 



+ 



f B fB*rn 2 B r gB*B^fn 2 B , + m|) 
4m b m B * ^{q 2 — m 2 B ,){q2 — m 2 B ) 
9b*Btt , 3/+ — /_ 



+ 



+ 



m 



B* 



it 



m 



B 



(2.8) 



where f B , /b*, f+ and /_ are defined by 



<0|A M (0)|S(p) >=i Pfi f B 
<0\V^0)\B*(e,p) >=^f B *m B , 
< n-(q)B°(q 2 )\B*-( qi ) >= (f + P, - f_q^ . 
Similar definitions hold for the D* — > Dn decay. 



(2.9) 

(2.10) 
(2.11) 



The second (II) integration region in (|2.7|) is defined as follows: ml < s < sq and 
s' > s or ml < s' < Sq and s > Sq. In this region we can distinguish two contributions. 
The first one is obtained by coupling the vector current to the pion and the B and the 
second one arises by coupling the pseudoscalar j'5 current to 7r and B*. The first term is 
as follows 

Aii{0,q l ,q 2 ) - 



2m b q\ - m\ 



3F 1 -mF -F 1 



(2.12) 

describing the coupling of the vector 



where F and Fi are the usual form factors 
current to the B and 7r. We can safely assume that F\ = F\{q\) is dominated by the 
B* pole ]19| , [20 , I^TJ , therefore the contribution from this form factor would be omitted 
to avoid double counting; assuming that the q\ dependence of F is given by the low lying 
+ pole (whose mass we denote by m) we obtain 



A n {^qlq 2 2/ 



f B m 2 B m 2 F (0) 



2mi, 



q\ — m 2 B q 2 — m 2 



(2.13) 



As for the second term in the II region, it would be proportional to the ^0(02) form factor 
(in the notations of ref. |18|). Assuming again the dominance of the low lying 0~ 5-pole 
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we would obtain a term which coincides with Aj and therefore it will be omitted to avoid 
double counting. 



Let us finally consider the third (III) region in the (s, s') plane, as defined by s, s' > s . 
Assuming duality, in this region Am should coincide with the asymptotic (qf 
— oo) limit of ( |2.3| ); therefore the spectral function is given by: 
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Pm(s,s') 



m b f n - 



< uu > 
3/. 



it 2 6(s - s )9(s' - s )S(s 



(2.14) 



Let us now discuss, before taking the single Borel transform of A(0,q 2 ,q%), the role of 
the terms (|2.13|) , (|2.14|) and the single pole contributions to Ai in (|2.8| ) (in the literature 
they are called "parasitic" or non-diagonal terms). In general these terms contain new 
unknown quantities as well as several uncertainties; however we can exploit their different 
dependence on the Borel parameter M 2 as compared to the resonance term (the first 
contribution in fl2.8| ) ); as a matter of fact, putting q\ = q\ and taking the Borel transform, 
one gets the following sum rule: 



7 

M 2 



+ A(M 2 ) = exp(fi/M 2 ) ( 



m b f n - 



< uu > 
3/. 
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+ 
+ 
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1 

M 2 
1 

2M 1 



2 < uu > ttl 2 lOm^m 2 m 2 ) < uu > 



3/* 
9 



3/. 



+ 



Wq < uu > ml 



6f, 



+ 



mQ < uu > ml 



f B f B * m\ gB*Bn 

-im^ms* 



(3m|» + m 2 B ) 



and Q = (m 2 B + m 2 B * 



2m 2 )/2, while A(M 2 ) is given by 
A(M 2 ) = d + d ie - 5/M2 



(2.15) 



(2.16) 



(2.17) 



with 

S = s -m 2 B (2.18) 
and do, d\ unknown quantities. We note that A(M 2 ) takes into account all the possible 



parasitic terms, either deriving from eqs. (|2.13| ),( p. 14[) or from the last two terms on the 
r.h.s. of (pTB"!). We also note that in deriving (|2.17|) and (|2.18 
the usually accepted values of s (s — 33 — 36 GeV 2 for the B 



we put m 2 



2 , 



So, which, for 
6-8 GeV 2 for the 



charm case) is a reasonable approximation. Moreover we put m B 



m B *)/2 and 



D 



(m 



D 



+ m 2 D *)/2. 



We can now differentiate (|2.15| ) in the variable 1/M 2 so that we remain with a new 



sum rule (NSR) for the variable 7 + d\/d(l/M ). Since 7 does not depend on M , we 
differentiate again this equation, so that we obtain an equation in the unknown quantity 
d\. Substituting this expression in NSR we finally obtain the following sum rule for g B * Bv : 
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mi 



M 2 



m^ < uu > 



mi 
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, 2 + m^ 



3ml 



lOf^mjml 
95 



+ 



m\ < uu > ml 



OfJ 



1 + 



6ml 
M 2 



(2.19) 



The analysis of the sum rule follows the usual criteria adopted in the literature: we check 
the existence of a region in M 2 where the OPE displays a convergent structure and where 
the highest states give a small contribution, as we discuss in the following section. Using 
mj, = 4.6 GeV and m c = 1.34 GeV |E| for the heavy quark masses we obtain, in the 
region M 2 > 30 GeV 2 for B and M 2 > 8 GeV 2 for D: 



/b Ib* 9B*Btt 

Id Id* 9d*Dtt 



0.56 ± 0.10 GeV 2 
0.34 ± 0.04 GeV 2 



(2.20) 
(2.21) 



the uncertainty is mainly due to the variation of the threshold Sq. It should be noticed 
that the parasitic terms represent a relevant fraction (40 — 50%) of the sum rule. 



The strong couplings gp*pK can be derived from ( j2.20| , |2.21| ) once the leptonic constants 
are known. Some authors (13], [14| suggest that in eqs. ( |2.15|) one should use as an input the 
leptonic constants obtained by two-point QCD sum rules at order a s = 0, to be consistent 
with the calculation of the correlator ( |2.2|) carried out at the same order in a s . Using 
f D = 170 ± 10 MeV, f D * = 220 ± 24 MeV, f B = 150 ± 20 MeV, f B * = 190 ± 10 MeU, 
we obtain: 



9B*B-it 
9d*Dtt 



20 ±4 
9 ± 1 . 



(2.22) 
(2.23) 



It should be noticed, however, that in a sum rule the 0(a s ) terms appear in the coefficients 
of the condensates, where they usually give small corrections, and in the perturbative 
contribution, where they are sizeable in the case of the leptonic constants. On the other 
hand, in the determination of gp*p 7T the perturbative contribution (given by the operator 
1 of the expansion of the correlator in eq. ( |2.2|) calculated between the vacuum and the 
pion state) vanishes at all orders in a s ; this is evident if one considers two-point functions 
in the external pion field 



13 



[14|. For this reason it seems to us that a s corrections 
could be included in the calculation of the leptonic constants, where they are relevant, 
whereas they could be safely neglected in all other quantities. Using the results obtained 
by two-point function QCD sum rules (with 0(a s )corrections) fp = 180 ± 30 MeV, 
f D = 195 ± 20 MeV Q and f B * = 213 ± 34 MeV, f D * = 258 ± 26 MeV obtained by a 



direct calculation in [^3| and by an analysis of the heavy quark spin symmetry breaking 
in p| we get for g P * P7T : 



9b*b-k = 15 ±4 
9D*Dn = 7 ±1 . 



(2.24) 
(2.25) 



For a review see 12, E2 
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Using the value in eq. flggSJ) we predict: T(D* + -> L>°7r+) = 10 ± 3 KeF. Moreover, 
assuming the dominance of the 5* pole in the form factor Fi(q 2 ) for B — > 7r£i/, we obtain, 
using ( p. 24 ): i*i(0) = 0.30 ± 0.07 in agreement with other QCD sum rules calculations of 



this form factor p], 0. 

Let us discuss an interesting consequence of our result. In a recent paper |23 it has 



been suggested that the scaling laws of the heavy quark effective theory can be used to 
derive the ratio \V u b/V c d\ from the comparison of the spectra B —>■ tt£u and D — > islv near 
the zero recoil point. A crucial ingredient in this strategy is the ratio 



Rbo = pR^l^L (2 . 26) 
V m B Id* 9d*Dtt 

which is equal to 1 in the infinite heavy quark mass limit. Our results point to a value: 

R BD = 1.06 - 1.11 (2.27) 

depending on the role of a s corrections (the highest value corresponds to a s = 0); this 
O(10%) correction should be taken into account when determining V uh in this approach. 



3 The limit ttiq — > oo 



The infinite heavy quark mass limit (mg ^ oo) in eq.( [2.15 ) can be performed according 
to the procedure already applied to leptonic constants [^], |f27| and form factors In 
terms of low energy parameters the quantities in eq.( [2.15|) can be written as follows: 



m B 

m B * - m B 



O 



1 

m b 



B 



B* 



(3-1) 



uj represents the binding energy of the meson, which is finite in the limit mj — > oo; 
the scaled leptonic constant F is independent of (if one neglects a small logarithmic 
dependence). F has been computed by QCD sum rules: for to = 0.625 GeV and the 

threshold y Q = in the range 1.1 - 1.4 GeV the result is F = 0.30 ± 0.05 GeV 3 / 2 (at 

the order a s = 0) p4] and F = 0.41 ± 0.04 GeV 3 / 2 (including radiative corrections) p7| . 



The sum rule for g in eq.(L5) can be derived from ( |2.15| ) after having expressed the 
Borel parameter M 2 in terms of the low energy parameter E: M 2 = 2m\,E. One readily 
obtains: 
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~3fT + 
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(3.2) 
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This sum rule must be studied in the region of the external parameter E where the 
OPE is assumed to converge and where the contribution of higher resonances is small 
("duality" region). By imposing the constraint that the various terms of the OPE in 
( |3.2| ) display a hierarchical structure, according to the dimension, we get that E should 
be larger than 1 — 2 GeV. As for the contribution of higher resonances, it is controlled 
in ( |3.2| ) by the parasitic term. It can be shown that this contribution, for values of the 
parameters in the range chosen above, is numerically of the same size of the contribution 
of the lowest lying states for any value of E. A prediction for g can be obtained by 
studying the stability plateau for (|3.2|), which corresponds to E > 4 GeV . 



We obtain: 

F 2 g = 0.035 ± 0.008GW 3 . (3.3) 
Therefore at the order a s = our result is: 

# = 0.39 ±0.16. (3.4) 

As we have discussed already, we could introduce a partial set of 0(a s ) corrections, i.e. 
those induced by F, which should represent the largest part of such corrections; in this 
case we get: 

g = 0.21 ±0.06. (3.5) 

The difference between (|3.4j) and ( |3.5| ) reflects the well known important role of radiative 
corrections in the determination of by QCD sum rules in the mg — > oo limit ||28|| . 
However it should be noticed that the value of g extracted from ( |1.5| ) using ( |2.22j ) or 
( |2.24| ) is in the range 0.19 — 0.25, which, if 1/rriQ corrections are not anomalously large, 



points to the smaller value in 



4 Conclusions 

The strong couplings gp*p n play an important role in heavy meson phenomenology. They 
are directly related to the decays of P* into P + ir and they are expected to be important 
for the semileptonic decays of P into a pion. They are also important inputs in the 
effective chiral lagrangians for heavy mesons. We have estimated the coupling constants 
9d*Dtt and gs*B-K in the soft pion limit and for finite heavy quark mass. Our calculation 
uses QCD sum rules as applied to the correlator of a vector and a pseudoscalar current 
between pion and vacuum. It is done for b and c states, thus allowing to examine the 
validity of the infinite heavy quark mass limit. We have considered the contributions of 
all the operators up to dimension five to the operator product expansion: it turns out that 
the higher dimension ones give small contributions. The soft pion limit and the single 
Borel transformation give rise to parasitic terms which are not exponentially suppressed. 
They give an essential correction to the value of the coupling, of the order of the lowest 
resonance contribution. The duality region, where the OPE shows convergence and the 
higher resonances give small contributions, is found for values of M 2 , the Borel parameter, 
higher than those usually found in the ordinary sum rules for semileptonic decays of heavy 
mesons. This could indicate a slow convergence of the OPE; nevertheless we have checked 
that in this region the parasitic terms do not give a contribution larger than 50 % to the 
sum rule. 

As we discussed before, the absence of the pure perturbative term in the OPE, where 
usually radiative corrections are relevant, seems to indicate that the largest a s corrections 
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in the sum rule could be taken into account by using the radiative corrected values of the 
leptonic decay constants. In the case of finite heavy quark mass, the difference between 
the two choices (i.e. leptonic decay constants with or without a s corrections) is within 
the uncertainties of our approach, while in the limit of infinite heavy quark mass there 
is a factor of two. Our numerical calculation of the sum rule, in the case when strong 
radiative corrections in the values of the leptonic constants are included, leads to values 
9b*b-k = 15 ± 4 and go*Dn = 7 ± 1. Without including these radiative corrections we are 
lead to gB*Bw = 20 ± 4 and Qd*d-k = 9 ± 1. The former case leads for the weak vector 
form factor between B and 7r, iq(0), to a value of 0.30 ± 0.07, in agreement with previous 
calculations. 
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